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ABSTRACT 

We present a Lagrangian description of the SU(2)/U(1) coset model per- 
turbed by its first thermal operator. This is the simplest perturbation that 
changes sign under Krammers-Wannier duality. The resulting theory, which 
is a 2-component generalization of the sine-Gordon model, is then taken in 
Minkowski space. For negative values of the coupling constant gr, it is classi- 
cally equivalent to the 0(4) non-linear a-model reduced in a certain frame. 
For g > 0, it describes the relativistic motion of vortices in a constant external 
field. Viewing the classical equations of motion as a zero curvature condition, 
we obtain recursive relations for the infinitely many conservation laws by 
the abelianization method of gauge connections. The higher spin currents are 
constructed entirely using an off-critical generalization of the Woo generators. 
We give a geometric interpretation to the corresponding charges in terms of 
embeddings. Applications to the chirally invariant U (2) Gross-Neveu model 
are also discussed. 
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1 Introduction 



Integrable perturbations of 2-dim conformal field theories (CFT) have been studied ex- 
tensively over the last few years. Zamolodchikov proved that there are perturbations of 
the form 

S = S C FT + 9 I $ (1.1) 



which take us away from the critical point, but the resulting theory still posesses an 
infinite number of local conservation laws [1, 2]. $ is typically a local spinless field with 
conformal dimension A$ which is found in the operator algebra of the unperturbed theory 
Sqft- The coupling constant g is a dimensionful parameter of weight (1 — A$, 1 — A$). 
In these cases, away from criticality, the chiral conservation laws of the CFT are replaced 
by 

dA s + dB s = (1.2) 

for appropriately chosen non-chiral currents {A s } and {B s }. Their existence is deeply 
related to the null vector conditions on the primary field $ driving the perturbation. 
Then, the integrals of motion off-criticality are 

Qs = J dzA s -dzB s , (1.3) 

with s ranging over an infinite set of values, depending on the model. 

In this paper we study the geometry and the classical conservation laws of the simplest 
Zjv-invariant CFT, ie the SU(2)/U(1) parafermion coset model, perturbed by its first 
thermal operator <3> = e±. In the large N limit, where all our work will be concentrated, 
this theory is described by the action 

s = — — + g(i-\u\ 2 ), (1.4) 

J 1— 7/ 



with | u \ 2 < 1. The first term is the usual classical action for the SU(2)/U{\) coset 
model and can be obtained from the SU{2) gauged WZW model after performing the 
neccessary gauge field integrations [3, 4, 5]. The potential term, which takes us away 
from the critical point while preserving the U(l) invariance of the theory, has conformal 
dimension (= limN^oo(2/N + 2)) and it corresponds to the first thermal operator of the 
unperturbed model. Many aspects of this model have been studied before, in connection 
with the relativistic theory of vortices in superfluids and other problems in field theory. 

The classical equations of motion that follow from the action (1.4) are 

ddu + UduBu 2 + g u (i _ | u |2) = o, (1.5) 
1 — | u I 

ududu ,o, / \ 

ddu + 2 +g u(l - | u | 2 ) = 0. (1.6) 

1— 71 
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This theory provides a 2-component generalization of the sine-Gordon model. Indeed, 
for u — u — cos9, the classical equations of motion reduce to 

2dd6 = g sin26 (1.7) 

and so the integrability of the full theory should generalize the results we already know 
for the sine-Gordon model. 

Integrable deformations of Zjv-symmetric models of CFT have been studied exten- 
sively [6, 7]. However, the Lagrangian description and the geometric interpretation of 
these perturbations have not been addressed in all generality. We focus on the sim- 
plest perturbation that changes sign under Krammers-Wannier duality and examine its 
geometry, in the context of perturbed parafermion models. The recent developments 
in the geometric interpretation of various CFT coset models as exact theories of black 
holes [5], provide the main motivation for considering this question with a perturbation 
switched on. We also write the classical equations of motion as a zero curvature con- 
dition with spectral parameter. Then, the abelianization method of gauge connections 
[8, 9] is employed to construct classically the infinitely many local conservation laws of 
the theory away from criticality, in a systematic way. Generalizations to other perturbed 
coset models are possible, but go beyond the scope of the present work. 

The Lagrangian approach to perturbed CFT coset models has been proven useful in 
other occassions for studying the integrability aspects of some special operators. For 
example, when a relevant perturbation by the (1,3) operator is applied to the minimal 
models of the Virasoro algebra with c < 1, the perturbed system is effectively described by 
the sine-Gordon equation (1.7), for appropriately chosen values of the coupling constant 
g [10]. In this case, the local conservation laws can be constructed systematically from the 
zero curvature formulation of the problem and they are related, as it is well known, with 
the Hamiltonian densities of the KdV hierarchy (ie, their flows are mutually commuting). 
More generally, affine Toda theories have provided a Lagrangian framework for looking 
at the problem of various integrable perturbations (see for instance [11]). 

The perturbed theory (1.4) has an interesting physical and geometric interpretation, 
when it is defined on Minkowski space. Based on some old work by Lund and Regge [12, 
13], we find that for g < its classical equations of motion describe a reduced form of 
the 0(4) non-linear a-model in two dimensions (see also [14]). For g > 0, it describes 
in a certain gauge the relativistic motion of vortices in constant external field. To put it 
differently, the physical picture for g > is that of a 4-dim bosonic string propagating 
in an axionic background of vortex type. At the conformal point g — 0, it describes 
the (transverse modes of the) free Nambu-Goto string in 4-dim Minkowski space in the 
orthonormal gauge (see also [15]). 

We should mention for completeness (and as independent motivation) that the per- 
turbed parafermion theory (1.4) is also interesting for the classical problem of two mass- 
less Fermi fields with contact 4-fermion interaction (ip a ip a ) 2 — (^J a 'j5^J a ) 2 [16]. The 
Nambu-Jona-Lasinio model in 2-dim Minkowski space (or Gross-Neveu model) provides 
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a chirally invariant generalization of the (multi-component) Thirring model. Exploiting 
the symmetries of this model, we may reduce it in a frame where 



WW = 9i , WW = 92, (1.8) 

waft - vwr = -2< gi 2 , ^r^2 - r^r = * (1.9) 

with (71 and 02 being constant. Here (1) and (2) denote the upper and lower components 
of the Lorentz spinors ip a , while summation over the fermion species a — 1, 2 is implicitly 
assumed. Then, it can be verified [17] that the composite complex fields 

u = -^j:r(i+ l5 )r , u = —±—Y,r{i- lb )r (1.10) 

satisfy eqs.(1.4), (1.5) with 

9 = 9\92- (1.11) 

If there is an arbitrary coupling constant in the 4-fermion interaction, it will enter mul- 
tiplicatively in (1.11). Consequently, in the absence of the 4-fermion interaction, the 
unperturbed SU{2)/U{\) coset model is recovered. 

This result is analogous to the well known relation between the Thirring and the sine- 
Gordon models [18]. Note in the present case that the chiral invariance of the Gross- 
Neveu model manifests as rotational (U(l) invariance) of the perturbed parafermion 
model. Indeed, in the variables (1.10), the chiral 4-fermion interaction is simply | u | 2 . 
Breaking chiral invariance of the theory by adding a mass term ijj a ifj a , corresponds to 
introducing the U (1) violating term u + u in the perturbed parafermion model. A useful 
consequence of our work is that the infinite many conservation laws of the theory (1.4) 
can be applied directly to the classical U(2) fermion system, using the transformation 
(1.10) in the special frame (1.8), (1.9). To the best of our knowledge, the systematic 
construction of the infinite many conserved charges of the Gross-Neveu model has not 
been carried out in detail. Although its integrable properties have been studied in general, 
the (off-critical) structure of its currents has not been recognized so far. 

Having presented an outline of the main issues and motivations of the present work, 
we describe briefly the organization of the remaining sections. In section 2, we examine 
the perturbation of the parafermion coset by its first thermal operator, in the framework 
of gauged WZW models. We also discuss the 1-soliton (and anti-soliton) solution that 
exists in the large N limit of the theory. In section 3, the physical interpretation and 
geometry of the classical theory (1.4) are subsequently described in Minkowski space 
for g > 0, g < and g = 0. In section 4, the infinitely many local conservation 
laws are obtained by the abelianization method of gauge connections. It turns out that 
for generic values of the coupling constant g, the higher spin currents of the theory 
are written in terms of an off-critical generalization of the generators. In section 5, 
connections with KdV type equations are presented. In particular, the conserved densities 
of the perturbed model are identified with the Hamiltonian densities of the (2-boson) 
KP hierarchy, thus generalizing the connection between the sine-Gordon equation and 
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SL(2) KdV hierarchy. Finally, in section 6, we present the conclusions and directions for 
future work. A geometrical interpretation of the conserved charges is also given in terms 
of embeddings. 

2 WZW Description of the Model 

We review first some standard results from the theory of Z N parafermions [19], in order 
to explain the form and the properties of the classical action (1.4). It is known that 
the field space of the SU(2) N WZW model (see also [20, 21]) contains N + 1 invariant 
fields with j = 0, 1/2, 1, 3/2, • • • , N/2; $ (0) coincides with the identity operator. 
Each is an SU(2) x SU(2) tensor with (2j + l) 2 components §m,m with to, to — 
—j, — j + 1, • • • , j — 1, j. These fields have conformal dimension j(j + 1)/ (N + 2) and 
to, to are the U(l) charges of §m,m m the two chiral sectors of the theory. The principal 
fields <j>2m,2rh °f the Z N parafermion theory are related with the components of the WZW 
invariant fields by 

where x an d X are the U(l) bosons of the two chiral sectors that are moded out in the 
construction of the SU(2) N /U(1) coset model. In this notation, <j k = <f>^J. are the spin 
variables and jj, k = 4> k k ^_ k are the dual spin variables of the parafermion theory. 

We are interested in the primary field with j — 1, in which case the 9 components of 
$W can be naturally identified [20] with the matrix elements 

^Trig-'Tagn), (2.2) 

where g, T are SU(2) group elements and Lie algebra generators respectively. The 
thermal operators (sometimes also called energy operators) are defined to be the U(l) 
neutral fields tj = (j)^ and so 

ei = Tr{g~ l a^ga^). (2.3) 

In the large iV limit, the conformal dimension of all e, goes to zero. An important 
property of the thermal operators under the Krammers-Wannier duality a <-> ji is 

e 3 - ( 1 ) J ej. (2.4) 

The Krammers-Wannier duality generalizes to arbitrary iV the relation between the low 
and high temperature phases of the Ising model. As we will see later, the change of sign 
in ei under a (/, implies, upon analytic continuation in Minkowski space, a duality 
relation between the 0(4) non-linear a-model and the relativistic theory of vortices in 
a certain frame. 
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Consider now the Lagrangian description of the SU(2)/U(1) coset model in terms of 
the SU(2) gauged WZW model [3, 4, 5]. Gauging the diagonal U(l) subgroup of SU(2) 
we obtain the action 

N f - 
S = S wzw + ^jTr (tAdgg- 1 - lAg^dg + AgAg' 1 - AA) , (2.5) 

where the gauge fields A, A take values in U(l). The classical equations of motion for A 
and A are in component form 

1 



2(1 - m 33 ; 

1 



Trig^dg a 3 ), (2.6) 



A = ~ 2(1^3) ^fe-Va), (2.7) 
where M 33 is given by 

M 33 = \ Trig- l a 3 gaA. (2.8) 
We fix the gauge by choosing SU(2) group elements of the form 

/ 9o + m igi \ 
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(2.9) 



V w\ go - igz J 



with g 2 + gi 2 + g 3 2 = 1. Then, in this unitary gauge, substituting the classical equations 
of motion for A and A in (2.5), we obtain the action of the SU(2)/U(1) coset model 

N r dudu + dudu , , 

S = to J l-\u\ 2 ( ° } 

in terms of the complex variables 

u = 9o + i93, u = g -ig 3 . (2.11) 

2 

Clearly, we have the condition | u \ < 1. 

The classical action (2.10) differs from the ordinary 0(3) non-linear a-model in that 
the target space metric is (1 — | n | 2 ) instead of (1 — | u | 2 ) . As a result, the topology 
of the target space is not that of a round sphere but of two bell touching together at 
the rim | u \ 2 = 1. The action (2.10) is valid only classically and therefore it provides 
a Lagrangian description of the parafermion theory in the large iV limit. Quantum 
mechanically there are 1/N corrections to the target space metric [22] and there is also 
a dilaton field from the path integral measure that insures conformal invariance at the 
critical point. In this paper the quantum mechanics of the problem will not be considered 
at all. 

The perturbation of Swzw by the first thermal operator e\ can be gauged similarly, 
with no extra effort, because e\ is the neutral component of the primary WZW field 
In the unitary gauge we have 

ei = 2 M 33 = 2 (2| u | 2 - 1) (2.12) 
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and so the action we obtain in the large N limit of the perturbed theory is essentially given 
by eq.(1.4), where g is the coupling constant of the perturbation. In deriving (1.4) we 
have shifted the action by a constant. This adjusts the zero of the energy density and has 
no effect on the classical equations of motion. We also note that unlike the sine-Gordon 
model (1.7), the coupling constant g of the 2-component generalization (1.5), (1.6) can 
not be made positive always. There are two phases in the perturbed theory, one for g > 
and one for g < 0, related to each other by the Krammers-Wannier duality (2.4). The 
conformal (critical) point g = is self-dual. Generalization of this formalism to other 
pertubations is also possible. It would be interesting to study the Lagrangian description 
and the geometric interpretation of perturbations driven by the higher thermal operators 
of parafermions, in a systematic way. The Lund-Regge formalism that will be adopted 
in the next section might be useful for handling the more general case as well. We hope 
to address these problems in the future. 

At the classical level, the coupling constant g can be normalized to 1 or —1, depending 
on its sign, with no loss of generality. This is possible because g is dimensionful of weight 
(1, 1). In Euclidean space, the two cases are related to each other by Krammers-Wannier 
duality, as it has already been pointed out. In Minkowski space, the sign of the coupling 
constant can change by interchanging the role of space and time coordinates. From 
now on we consider the theory (1.4) defined in Minkowski space and study its physical 
interpretation for g > 0, g < and g = in the light-cone coordinates 

d = d a + d T , d = d a -d T . (2.13) 

The space and time variables will be a and r in all cases. 

An important property of the perturbed parafermion theory is that it admits soliton 
solutions, in analogy with the sine-Gordon model. To make the comparison easier, we 
introduce the variables A and 9 

u = cos9 e iX , u = cosO e- iX . (2.14) 

In the static limit d = 3, the 1-soliton solution of the classical equations of motion (1.5), 
(1.6) was found long time ago [12]. For g > (and normalized to 1) the 1-soliton is 



9(z) = sin^ 1 I ^ A2 . - 1 = sin' 1 {Vl - A 2 sin (2 tan' 1 exp(zy/ 1 - A 2 )) } , 

I cosh iz v^l J^P 1 1 I 

(2.15) 

X(z) = Aj Z dz' tan 2 9(z'), (2.16) 

where A is constant. In the limit A — > we have A = and the usual sine-Gordon 
1-soliton solution 9 = 2 tan~ 1 expz is obtained (see for instance [23]). The corresponding 
anti-soliton solution is obtained by shifting 9 — > tt — 9. As we will see later, the parameter 
A is the U(l) charge of the solution and so the sine-Gordon model describes the zero 
charge sector of the theory. A localized lump travelling with constant velocity is obtained, 
as usual, by Lorentz transformation. 
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Multisolitons solutions can also be constructed, using standard techniques from soliton 
theory. The difference between g > and g < in the soliton solutions has been 
considered in ref. [24], together with some related issues. An analogous solution, for 
(7 = 0, was found by Bardacki et.al.[4]. In the static limit of the theory, 



and X(z) also given by eq.(2.16), solve the classical equations of motion without the 
thermal perturbation. 

The existence of soliton solutions in the perturbed parafermion model might have 
important consequences for the associated scattering problem. This sector has not been 
included so far in determining the S"-matrix of the theoryQ Naive consideration of the 
large N limit shows that S = 1 for this model [6, 25]. It is not clear at this point whether 
the soliton solutions persist for finite N or whether they are characteristic of the large N 
limit. It might be possible to address this question by considering the 1/N corrections 
to the classical action (1.4) from the gauged WZW point of view. In any event, the 
scattering of these solitons should be investigated separately in the future. 

3 Geometry of the Thermal Perturbation 

In this section we describe geometrically the classical equations of motion of the perturbed 
parafermion model. The main subject of the section is the embedding of a 2-dim surface 
S with local coordinates a, r and metric 



in a 3-dim space of constant curvature, which in turn is embedded in 4-dim flat space. 
The differential equations that the metricf] g and the extrinsic curvature tensor K of S 
have to satisfy, in order to have a solution to the embedding problem, are given by the 
Gauss-Codazzi integrability conditions (see for instance [26]). These conditions admit 
a field theoretic interpretation which helps us to understand the physics and geometry 
of the first thermal perturbation of the parafermion coset. The technical details of the 
embedding will be considered later and they are intimately related with the physical 
interpretation of the perturbed coset model. The cases g < and g > are treated 
separately, since the results we obtain depend crucially on the sign of the coupling con- 
stant. In this framework we also revisit the critical point g = and find that the 
SU(2)/U(1) coset model provides an effective description of the transverse modes of the 
4-dim Nambu-Goto string in the orthonormal gauge. 

The formalism we adopt here originates in the work by Lund and Regge [12, 13] (see 
also [14]), but since it is rather unknown we review briefly the main ideas applicable to 

*I thank T. Hollowood for some discussions on this point. 

t Although we use the same symbol g for the metric and the coupling constant, the distinction between 
the two will be obvious from the context they appear. 




(2.17) 



ds 2 = cos 2 8da 2 + sin 2 9dT 2 



(3.1) 
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our model. The considerations are entirely local, with no proper reference to boundary 
conditions. Throughout this section, the perturbed parafermion model is defined on 2- 
dim Minkowski space. Analytic continuation to Minkowski space is neccessary, in order 
to describe the physics and geometry of the theory via embeddings. 

(i) g < 0: In this case, the perturbed model (1.4) is classically equivalent to the 2- 
dim 0(4) non-linear a-model, reduced in a certain frame. Recall that the 0(4) model 
consists of four scalar fields r)i(a, r) (i = 1, 2, 3, 4) interacting through the quadratic 
constraint 

( Vl ) 2 = 1. (3.2) 
The solutions of the classical equations of motion 

(d a 2 - d T 2 ) Vt + ((d aVk ) 2 - (d rVk ) 2 ) Vl = 0, (3.3) 

describe 2-dim surfaces (parametrized by o and r) embedded in the 3-dim sphere (3.2), 
which in turn is embedded in flat 4-dim Euclidean space. The embedding variables are 
simply rji. It is natural to reduce the 0(4) non-linear a-model in the frame 

(£U) 2 + (d rVt ) 2 = 1 , (cU)(cU) = (3.4) 

by exploiting the invariances of the classical theory. When the conditions (3.4) are 
satisfied, the induced metric on the 2-dim embedded surface is given by eq.(3.1), where 

cos 2 6 = (d a r]i) 2 , sin 2 6 = (<9 T ^) 2 . (3.5) 

Note that the 0(4) model is defined in Minkowski space, while the associated embedding 
problem is entirely Euclidean. 

The formulation of an embedding problem requires apart from the metric g, knowledge 
of the extrinsic curvature tensor K. Its components are defined to be 

K aa = (d a \)Z^ , K TT = (d T \)Z^ , K aT = K TU = (d^d^zP, (3.6) 

where = Z^ x with 

We also choose Z^ = rji to complete an orthonormal tetrad in 4-dim Euclidean space, 
where the embedding takes place. The classical equations of motion of the reduced 0(4) 
non-linear a-model are equivalent in the special frame (3.4) to 

K aa = K TT . (3.8) 

This is obtained by multiplying the classical equations of motion (3.3) with Z^ 3 \ Mul- 
tiplication with Z^\ zf^ and z\^ leads to identities in the special frame (3.4). The 
reduction of the 0(4) model in this frame is analogous to the reduction of the chirally 
invariant U (2) Gross-Neveu model, described in the introduction. 
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It is clear that the classical equations of motion of the reduced theory can be formu- 
lated as an embedding problem in the 3-dim sphere (3.2), with the extrinsic curvature 
satisfying the physical requirement (3.8). This embedding is possible provided that the 
Gauss-Codazzi conditions on g and K are satisfied. In the present case they assume the 
form 

d T (tan6 K aa ) = d a {tan6 K UT ) , d a (cot6 K au ) = d T {cot9 K ar ), (3.9) 
R = (g al3 K a/3 f - g^g^K^K.s + 2, (3.10) 
where R is the scalar curvature of the 2-dim metric (3.1) 

R= . I - (O-O)- (3.11) 

smO cos 9 

The constant term in eq.(3.10) is the curvature contribution of the 3-dim sphere (3.2) 
where the embedding takes place. The first condition in eq.(3.9) implies the existence 
of a scalar field A such that K aa = cotO d a \ and K aT = cot6 d T \. Introducing complex 
variables u,u as in eq.(2.13) and light-cone coordinates as in eq.(2.16), we obtain 

udu — udu udu — udu 
K aa + K UT = i , K aa - K aT = i . (3.12) 

2 I u I \ 1 — I u I 2 I u I \ 1 — I u I 



Then, it is straightforward to verify that the remaining Gauss-Codazzi conditions become 
the classical equations of motion of the perturbed parafermion model in Minkowski space, 
with g = —1. 

The reconstruction of S (and hence rji) can be done by solving the Gauss- Weingarten 
equations for the orthonormal vectors Z^ 2 \ and Z^\ These are first order linear 
differential equations, forming a system, whose compatibility conditions are provided by 
the Gauss-Codazzi equations (for details see [26]). The classical equivalence between the 
two theories, both defined in 2-dim Minkowski space, has been established with the aid 
of a Euclidean embedding problem. In this case, the curvature of the 3-dim sphere (3.2) 
determines the coupling constant (it turns out to be —2g) and hence g is negative. 

The physical interpetation of the perturbed theory with g > will be addressed 
differently, without interchanging a <-> r. Following Lund and Regge [12], we find that 
the phase of the theory with g > describes (in a certain gauge) the relativistic motion 
of vortices in a constant external field. We also adopt this picture in order to show that 
the transverse modes of the 4-dim Nambu-Goto string in the orthonormal gauge are 
effectively described by the conformal limit of the SU{2) /U{\) coset model. 

(ii) g > : Consider first the string action in 4-dim Minkowski space 

S = -N J ^-detg dadr + f J A^(X)d a X»d T X h 'dadr --J F^d 4 y, (3.13) 

where a, r are local coordinates on the string world-sheet E and detg is the determinant 
of the 2-dim metric on S. The first term is the usual Nambu-Goto action, while the 
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second term represents a self-interaction of Kalb-Ramond type [27]. A^ U (X) is an anti- 
symmetric tensor (axionic background) and F^ is defined to be 

F» = \ t^ xp d p A uX . (3.14) 

X M ((j, r) (/j, = 0, 1, 2, 3) are the embedding variables of the string in 4-dim Minkowski 
space with signature — h ++. When the spatial components of A pu are linear in X i 
(i = 1, 2, 3), ie 

Aij = eij k X k , (3.15) 

this action describes the relativistic motion of vortices (strings) in a superfluid [12] (see 
also [28] for the non-relativistic limit). In this case, F l is identified with the velocity of 
the fluid 

F i = e^djAko = v l (3.16) 

and the last term in eq.(3.13) represents the hydrodynamic action of the superfluid. 

Following Lund and Regge we study the relativistic motion of vortices in a uniform 
static external field. We choose a Lorentz frame in which 

X° = t , F i = (3.17) 

and F° is constant and introduce the effective coupling constant of the theory 

fF° 

°=w- (3 - 18) 

At this point we do not normalize c to 1, keeping its dependence on / and N explicitly. 
We need this in order to understand the physical interpretation of the limit c — ■> that 
will be considered later. We also choose the orthonormal gauge for the induced metric 
on the vortex (string) world-sheet, which leads to the quadratic constraints 

(d a X t ) 2 + {d T X % f = 1 , (d a X t )(d T X t ) = (3.19) 

in the Lorentz frame (3.17). This gauge choice is analogous to the reduction of the 0(4) 
non-linear cr-model (3.4) introduced earlier. Then, the vortex dynamics is determined 
entirely by the classical equations of motion 



(d a 2 - d/)Xi + 2 c ((d a X) x (d T X)) . = (3.20) 

plus the constraints (3.19). 

Note that for uniform static external field, the last term in the action (3.13) can be 
dropped out (it is just a number). The theory in this case describes a 4-dim bosonic 
string in an axionic background of vortex type. For / = 0, ie c = 0, it describes the 
propagation of a free bosonic string in 4-dim Minkowski space in the orthonormal gauge 
with Xo = t. It is important to realize that the Lund-Regge model for vortex dynamics in 
a uniform static external field is classically equivalent to the perturbed parafermion coset 
(1.4) with g = c 2 > 0, upon analytic continuation in Minkowski space. This formalism is 
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certainly unphysical for g < 0, since the coupling constant of the vortex self-interaction 
term will be imaginary in that case. We also point out that the results presented in the 
sequel are independent of the Lorentz frame X = r, up to Lorentz transformations in 
the (a, t) space. 

To illustrate the vortex-like description of the perturbed parafermion model with 
g > it is convenient to incorporate the classical equations of motion (3.20) in the 
integrability conditions of an embedding problem, in analogy with the g < case. Let S 
be the projection of the string world-sheet S in the X = r hyperplaneQ. S is a Euclidean 
surface, which in the orthonormal frame has an induced metric given by eq.(3.1). We 
have 

(c^X,) 2 = cos 2 6 , (d T Xi) 2 = sm 2 e. (3.21) 
In analogy with the g < case, we also consider the extrinsic curvature tensor of S, 

K aa = (d^X^Z® , K TT = (<9 r 2 X)^ {3) , K UT = K Ta = {d ff d T Xi)Z^ (3.22) 

with Z^ = Z« x Z®, where 

In terms of these variables, the classical evolution of vortices in the orthonormal frame 
is entirely determined by the embedding of S in the 3-dim Euclidean space Xq = r, 
which in turn is embedded in 4-dim Minkowski space. The components of the extrinsic 
curvature tensor have to satisfy the condition 

K aa - K TT + 2c sinOcosO = 0, (3.24) 

which follows from the classical equations of motion (3.20) by multiplication with Z\ . 
Multiplication with Z^ and Z^> leads to identities in the special frame (3.19). Hence, 
eq.(3.24) encodes all the information contained in eq.(3.20). 

The embedding problem here is different in that the 3-dim space where the embedding 
takes place is flat and the extrinsic curvature satisfies the condition (3.24) rather than 
(3.8). The Gauss-Codazzi integrability conditions have to be satisfied, however, in order 
to have a solution to the problem. In this case they read 

d T (tan9(K aa + K TT )) = 2d a (tan6K aT ) , d a (cot6(K a(T + K TT )) = 2d T (cot9K ar ), (3.25) 

R = (g al3 K aP f - sTg^K^K^ (3.26) 

where R is the curvature of the 2-dim metric (3.1), given again by eq.(3.11). Introducing 
light-cone variables z, z and complex coordinates u, u as before, we find that the first 
condition in eq.(3.25) implies 

udu — udu udu — udu 
K a(T +K TT +2K aT = i 2 , K aa +K TT -2 K aT = i = . (3.27) 

I U I \ 1 — I U I I U I \ 1 — I u I 



■^To know X, it is sufficient to consider S for X° = t. 
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It is staightforward to verify using eq.(3.24) that the remaining Gauss-Codazzi integra- 
bility conditions become the classical equations of motion of the perturbed parafermion 
model in Minkowski space, with g = c 2 . The reconstruction of S (and hence Xj) can 
be done by solving the corresponding Gauss- Weingarten equations for the orthonormal 
vectors Z« and . 

(iii) g = 0: When the coupling constant / of the antisymmetric tensor is zero, the 
theory reduces to the Nambu-Goto string propagating in 4-dim Minkowski space. As it 
has been pointed out already, in the orthonormal gauge with X = r, the unperturbed 
SU(2)/U(1) coset model, when it is defined in Minkowski space, describes the classical 
dynamics of the transverse modes of a free string. This result is very intriguing and 
probably it can be generalized to more arbitrary backgrounds. In a separate publication 
[15], the implications of this equivalence are analyzed in detail. The infinitely many con- 
servation laws of the parafermion coset model can be easily applied to string dynamics to 
yield hidden symmetries (Backlund transformations) in their classical solution space. In 
the next section we analyze the structure of the conserved charges for arbitrary values of 
the coupling constant g and establish their relation with Woo. Then, a geometric inter- 
pretation of the higher spin currents follows immediately from the embedding problem 
that describes the classical equations of motion of the theory. 

The classical equivalence between 4-dim string theory and the perturbed parafermion 
model with g > is not restricted to the special frame X = r. For X = f , where 

f = sinha a + cosha r , a = cosha a + sinha r (3.28) 

are related by Lorentz transformation in (a, r) space, the embedding problem remains 
essentially the same. The projection of the string world-sheet is performed now in the 
Xq — f 3-space, but the classical equations of motion that follow from the Gauss- 
Codazzi conditions remain unchanged [12]. For this reason, the equivalence between the 
two theories is independent of the Lorentz frame. This issue does not arise for g < 0, 
because the physics of the problem does not require the choice of a Lorentz frame to 
formulate the corresponding embedding. 

The soliton solution (2.14)-(2.15) of the classical equations of motion with g > 
becomes relevant for string propagation in an axionic background of vortex type. It 
would be interesting to study further the localized properties of this soliton solution, 
directly in the string variables Xj. The issue of boundary conditions and the physical 
interpretation of coset model solutions has to be addressed properly in the future, in the 
context of string theory. 

4 Local Conservation Laws 

We turn now to the explicit construction of the infinitely many conservation laws of the 
perturbed parafermion theory for arbitrary values of the coupling constant g. For this, it 
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is convenient to rewrite the classical equations of motion (1.5), (1.6) as a zero curvature 
condition. We introduce the linear system of differential equations 

<9$ = A $ , <9$ = B $, (4.1) 

where A and B are 2x2 matrices depending on z, z and a spectral parameter A 

A = A + \E , B = B + \- 1 B 1 . (4.2) 

The spectral parameter should not be confused with the scalar field A appearing in the 
definition (2.13) of u and u. A , B and B\ are taken to be 

\ 



A, 



1 



u 



'l-\u 



Bn = 



' (2| u | 2 - l){udu-udu) U\u \ 
^ 4i | u | -y/ 1 — I - 



u 



udu 



udu — (2| u | — — udu) j 



udu — udu 



9 
4 



u 
( 

\ 



E, 



2\u\-l 



—2i I u 



u 



|2 \ 



2i I -u 



(4.3) 
(4.4) 

(4.5) 



■(2\u\ 2 -l) ) 



and E is the constant matrix 



/ 1 



E 



\ 



(4.6) 

\0 -1 ) 

Then, it may be easily verified that the compatibitity (zero curvature) condition of the 
linear system (4.1), 

[d - A , d-B] = 0, (4.7) 

is equivalent to the classical equations of motion (1.5), (1.6), for all values of A. 

The physical meaning of the spectral parameter A can be understood in terms of the 
Lorentz transformation (3.28). Suppose that the linear system (4.1) is given initially 
with A = 1. The Lorentz transformation (a, r) — > (<r,f) amounts to the rescaling d — > 
exp(-a) d, d — > expa d and A — > exp(-a) A , B — > expa B . Absorbing the Lorentz 
factor, is equivalent to choosing A = expa in the gauge connections (4.2). Consequently, 
the independence of the classical equations of motion on A is really a statement about 
Lorentz invariance. As we will see shortly, the presence of a spectral parameter in the 
linear system (4.1) is neccessary to derive the infinitely many conservation laws of the 
theory in a systematic way. The results are valid equally well in Euclidean and Minkowski 
space. 

Consider the static limit d — d first. It follows immediately from [d — A , d — B] — 
that the quantities 



I n = \ Tr(A - BY 



n=l, 2, 3, 



(4.8) 
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are conserved, ie dl n = 0. This is a well known result for 1-dim integrable systems. If 
the number of degrees of freedom is finite, the quantities I n will not be all independent. 
For the present model we find that only I 2 is independent, because 

/ 2fc -i=0, I 2k = (I 2 )\ k = l, 2, 3,.-. . (4.9) 

Explicit calculation shows that 

/ 2 = ^ 2 -#-AQ + ^Q + ^ + A 2 + T !^, (4.10) 



where 



dudu l2 „ udu — udu . „ 

H=- ~-2+9 \uf, Q=- — 2 -. (4.11) 

1 — | u ] 1 — I u I 

Since dl 2 = for all values of A, we arrive at the conservation laws 

OH = = dQ (4.12) 

in the static limit of the theory. H is the energy and Q is the £7(1) charge of a given 
configuration u, u. For the 1-soliton (2.14), (2.15), with g normalized to 1, we find that 
H = 1 and Q = 2A, provided that the charge of the theory is made real by multiplication 
with %. This gives a physical interpretation to the arbitrary parameter A that determines 
the soliton solution. Solitons with zero U{1) charge are simply sine-Gordon solitons. 

In the general case, the conserved quantities of the theory are not given by eq.(4.8) 
and a field theoretical prescription is required to obtained the currents associated with the 
2-dim zero curvature condition (4.7). The abelianization method of gauge connections [8, 
9] provides the algorithm for their systematic construction. The main idea is to introduce 
a family of A-dependent gauge transformations 

CO 

T=l + X>-^, (4.13) 

i=i 

where ti depend on z, z and take values in GL(2) with det T ^ 0. Then, it is always 
possible to choose {ti} appropriately so that the gauge transformed connections 

CO 

T-\d - A - \E) T = d + A^ A, - XE, (4.14) 

3=0 

CO 

T~\d -B - X^Bi) T = d-J2 X~ j Bj (4.15) 

commute with the constant matrix E (4.6), ie 

[A 3 , E] —0 — [B 3 , E] (4.16) 

for all values of j. Since E is not the identity matrix, {A,} and {Bj} commute among 
themselves and the zero curvature condition (4.7) yields upon gauge transformation an 
infinite number of functionally independent non-trivial conservation laws 

BAj + dBj = , j = 0, 1, 2, • • • , (4.17) 
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one to each order in the 1/A expansion. The presence of a spectral parameter in the zero 
curvature formulation of the problem is indeed essential. The validity of the equations 
for all values of A implies order by order the infinite hierarchy of conservation laws 
(4.17). For E diagonal, as it is the case here, the abelianization method amounts to the 
diagonalization of the gauge connections A and B. It is important to emphasize that the 
diagonalization of B is not automatic and it can be achieved only on-shell. 

Having presented the essential ingredients of this method, we may proceed with ex- 
plicit calculations. Note first that the diagonalization of A implies an infinite set of 
conditions. We have 

A + [E , t 1 ]+A = (4.18) 
to order A , while for j > we obtain the recursive relations 

Aj + [E , t j+1 ] + A tj - dtj + ]T tj- k A k = 0, (4.19) 

fc=0 

one to each order in \~K Since all Aj have to commute with E, their most general form 
is 



A 3 



(hf \ 



V o hf J 



J > 0. (4.20) 



3 

The problem now is to solve these recursive relations, in order to determine the form of 
Aj and the gauge transformation T. 

Note that the solution to the diagonalization problem of A is not uniquely determined. 
There is the freedom to set the diagonal elements of all ti equal to zero, with no loss of 
generality. Indeed, any 2x2 matrix can be written in the form X + [E, Y], where X 
is a diagonal matrix and Y off-diagonal. Taking this decomposition into account in the 
Aj + [E, tj+i] part of the recursive relations (4.18), (4.19), it can be easily seen that for 
diagonal Aj, as they should be, t, + i can be chosen so that their diagonal elements are 
zero. This gauge choice fixes the form of U and Aj uniquely and the solution that results 
from the recursive relations has no free parameters; everything is functionally dependent 
on the matrix elements of A . This gauge choice will be assumed from now on. 

Gauge transformations with non-zero diagonal elements, modify the solution for Aj 
by total 9-derivative terms, via the recursive relations (4.19). This modification is sup- 
plemented in Bj by substracting the (^-derivative of these additional terms. This has no 
effect on the local conservation laws (4.17) and it reflects the freedom we have in writing 
down the corresponding currents. We also note that the gauge transformation T is not 
neccessarily restricted to SL(2), but it can take values in GL(2). In the gauge where the 
diagonal elements of all £j are set equal to zero, the determinant of T is 

detT = 1 - tftf\- 2 - {tftf + tjHl 2 )\- 3 + ■■■ . (4.21) 

Therefore, it is natural to expect that for j > 2 the matrices Aj will not be traceless. 
We could have chosen a gauge (order by order in 1/A) so that detT — 1, to preserve the 
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traceless condition on all gauge connections. However, the two gauges are related to each 
other by the U(l) element of GL(2), which is proportional to the 2x2 unit matrix, 

1 + I tftf\~ 2 + X - (tlHf + tf tl 2 )X~ 3 + • • • (4.22) 

and the difference on the diagonalized gauge connections A, B is therefore a total deriva- 
tive trace term. We choose to work with a gauge in GL(2), rather than SL(2), in order 
to simplify the form of the recursive relations. 

From eq.(4.18) we have immediately 

} = -K 1 = -4 2) (4-23) 

and 

t\ 2 = ~- 2 K 2 , t? = Ul\ (4.24) 

in terms of the matrix elements of A . This is the initial data for the recursive relations 
with j > 0. To bring the gauge connection A into the desired form, the following system 
of equations has to be iterated: 

hf = -AlHf , hf = -Aft?, (4.25) 
2 tf +1 = -AlHf + dtf - £ 4\ 12 -*, (4-26) 

k=0 

2 tf +1 = -AlHf - dtf + £ h?tf_ k . (4.27) 

k=0 

These relations follow from eq.(4.19), using the diagonal form (4.20) for Aj and the 
off-diagonal gauge for tj. 

The solution (4.23) for the matrix elements of A can be written in terms of u and u 

as 

i 1 2 

4 1} + 4 2) = 0, " -1 (udu-udu). (4.28) 

2 I u I (1 — | u I ) 

On the other hand, since _B = -Boj w e derive to this order the non-chiral conservation 
law 

dJ + dJ = 0, (4.29) 



where 

udu — udu 
1 — \u\ 



= 2(ft< 1) - ftf) - 3 (log f) , (4^30) 



J^ = «-fl + 9(l»» ")• ("I) 
1 — I u I \ UJ 

The A conserved currents are written conveniently in this form, in order to identify J 
and J with the two components of the U{1) current associated with the symmetry of 
the action (1.4) under u — > ue' lt and u — > -ue~* e . Note that the conservation law (4.29) is 
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2 

independent of the coupling constant g, simply because the term | u | in the potential of 
the action is invariant under U(l) rotations. In the static limit d — d, we have J = J = Q 
and the conservation law dQ = is recovered. 

Using eqs.(4.24) and (4.25), we obtain the following expression for A 1 in terms of u 
and u: 

h? + = , h? - hf = (4.32) 

1 — I u | 

Explicit calculation also shows that 

B 1 = B 1 + [B , ti] -dh = ^g (2| u | 2 - 1) £, (4.33) 

where the second part of the equation is valid only on-shell. Therefore, to order A -1 , the 
conservation law for the zz and zz components of the stress energy-tensor of the theory 
is obtained, 

*(iqvi I )+»«i «!' = »■ <« 4 > 

At the conformal point g — 0, this reduces to the chiral conservation law for the zz 
component of the stress-energy tensor of the SU(2)/U(1) coset model. In the static 
limit, the result dH = is also recovered. 

Iteration of the recursive relations for j > 2 yields higher order non-chiral conserva- 
tion laws, all depending on the coupling constant g. In the static limit they all reduce to 
the conservation of H and Q, but in the field theory case we are considering now they 
turn out to be functionally independent and hence new. To understand their nature, 
in general, it is convenient to compare them with the infinitely many symmetries of the 
SU(2)/U(1) coset model with g = 0. In the latter case, it is known that there is a chiral 
Woo symmetry, whose generators are bilinear in the parafermions of the model [29-31]. 
In particular, introducing the parafermion currents 



, , = ° v „ V + . , = = V ; 1.35) 

u 1 



where 

V± = exp (±- j dz J-dz J^j (4.36) 

are defined in terms of the non-chiral U(l) current of the theory, the W^o generators (in 
a quasi-primary basis and up to an overall normalization) are 



W s = V , , J , , m + d s - k - 2 ^ (4.37) 



with s = 2, 3, 4, • • • . For g — 0, ip± and W s are all chirally conserved, provided that 
the classical equations of motion of the SU(2)/U(l) coset model are satisfied. For g ^ 
this is not true, but it still makes sense to define parafermion currents ip± as in (4.35) 
and (4.36), because the potential | u | 2 is [/(l)-invariant. 
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The conservation laws that result from the diagonalization of the gauge connections of 
the perturbed SU(2)/U(1) coset model can be written systematically for all g, using the 
off-critical generalization of the generators (4.37) in terms of the parafermions ip±. 
Note that for g ^ 0, dW s can not be always brought in the form dX for appropriately 
chosen X . Although this is true for s = 2, for higher spin currents it is not so. To 
describe the conserved currents of the theory in terms of {W s } off-criticality, we have 
to introduce appropriate polynomial combinations of the generators, depending on s. It 
is the abelianization method that provides the algorithm for writing down these higher 
(non-chiral) conservation laws of the theory in the form (4.17). 

We find that for j > 2, the trace of Aj is a total derivative of currents, which are 
composed from the subleading componets hf^ — hfl^ hfl 2 — hfl 2 , etc. Hence, the only 



functionally independent conservation laws are obtained by considering hj — hj . For 
example, for the first few values of j (apart from the ones already discussed) we find 

h<P+W = -\d(hP-h?), (4.38) 

hV + h? = -\d{h?-h?\ (4.39) 
h? + hf = -1 d (24(/4 1} - hf) - (h? - h?) 2 - d 2 (h? - h?)) (4.40) 



and so on. As for the h^p — h!p components of Aj, the results of the calculation are 
summarized in the appendix for the first few values of j. The complexity of the expres- 
sions increases considerably for higher values of j, but they are all calculable order by 
order. Unfortunately, no closed expression for arbitrary j is available at the moment. 
The form of h!p — h^ indeed simplifies considerably, when they are written in terms 
of the parafermions off-criticality. Extending the definition of W s to all values of the 
coupling constant g and introducing the explicit dependence of A on u and u, we obtain 



after some lengthy calculation 

hP - h { ? = W 2 , (4.41) 

h? - hf = -1 W 3 , (4.42) 

- = 1 (W 4 + 5(W 2 ) 2 + ? d 2 W 2 ^ , (4.43) 

hP - hf = --i- (W 5 + 21W 2 W 3 + 6 d 2 W 3 ) (4.44) 



and so on. This summary shows that hj — hj are functionally related to the W s 
currents, as advertized. Note that these polynomial combinations are independent of g. 

The form of the infinitely many conservation laws of the theory (with arbitrary g) 
is determined completely, by performing the gauge transformation T on the gauge con- 
nection B as well. The iteration of the recursive relations also determines {t{\ order by 
order. Explicit results for the first few values of j are given in the appendix, together 
with Bj. It can be verified that the matrices Bj are all diagonal on-shell, as required 
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by the abelianization procedure. The final result for the local conservation laws of the 
theory (discarding total derivative terms and overall normalization factors) is 

dW 3 + g d{udu - udu) = 0, (4.45) 
d(W 4 + 5W 2 2 ) +g (d 2 | u | 2 + 5 ^ U j ~} duduj = 0, (4.46) 



d(W 5 + 21W 2 W 3 ) + gd\ ud 3 u - ud 3 u l -^-(dud 2 u - dud 2 u) 

V i ~~ I u I 

O Q I I \ 

+ 7 1 U I 2 {udu - udu)dudu = 0. (4.47) 



u 



r J 



The higher order conservations laws can be constructed in a similar fashion, but it is 
computationally difficult to find closed expressions for them, in general. Of course, we 
also have a complementary set of non-chirally conserved currents which is obtained by 
interchanging u <-> u and d <-> d. 

In the conformal limit g — > 0, the chiral W 7 ^ algebra is recovered, but in a polynomial 
basis. Its commutation relations are not linear in this basis, but we know that there exists 
the quasi-primary basis (4.37) where the linear structure of the algebra is manifest. For 
g ^ it does not make sense to talk about the algebra of the non-chiral currents, 
because the symmetries they generate are global; it is only at the conformal point that 
global symmetries are promoted to local. What makes sense, however, is to consider the 
corresponding charges for all s = 1, 2, 3, • • • , which are conserved and in involution. 
Since they are infinitely many of them, we have a rigorous method (and an algorithm) 
for establishing the complete integrability of the parafermion model perturbed by its first 
thermal operator, in the large N limit. 

The path ordered exponential 

P exp ^ 2 dzA + dzB^j (4.48) 

is independent of the path joining two space-time points P\ and P 2 , provided that A, 
B satisfy the zero curvature condition (4.7). Since the expression (4.48) depends only 
on the end points Pi and P 2 , we may use it to obtain an alternative description of the 
conserved charges, by considering paths joining two spatially far apart points at different 
times. The transformation to the diagonal gauge connections A, B makes the path 
ordering uneccessary. Then, the 1/A expansion of the new variables yields the infinite 
set of conservation laws we have already discussed. The generating function for the 
conservation laws can be obtained from the transition matrix of the theory, using the 
inverse scattering method (see [32] and references therein). This method has been applied 
to the model (1.4) by Lund [13] and Kulish [33], in order to construct the corresponding 
action-angle variables. In this regard, our results are complementary to theirs, but more 
explicit. The main point of this section was the realization that and the off-critical 
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generalization of its generators in the parafermionic representation (4.37) determine the 
structure of the higher spin conservation laws completely. 

The systematic construction of the conserved currents was based on the Langrangian 
description of the model and the zero curvature formulation of its classical equations of 
motion, rather than the null-vector conditions on the parafermion first thermal operator. 
The quantum mechanical generalization of these results should be straightforward, but 
computationally difficult. The quantum inverse scattering method could be used to 
understand the relation with the null-vector conditions in a systematic way. We think 
that the cohomological framework of Feigin and Frenkel [34] might be also appropriate for 
investigating further the quantum mechanical formulation of the abelianization method. 



5 KP-like Structure of the Currents 

The conservation laws of the perturbed parafermion theory can be described systemati- 
cally in terms of the KP hierarchy. We claim that the currents h!p — \ with j > 1, can 
be interpreted as Hamiltonian densities of the KP hierarchy, thus providing a more direct 
way for their explicit construction. This generalizes the well known relation between the 
conserved densities of the sine-Gordon model and the SL(2) KdV hierarchy, to the full 
coset theory. 

Recall that the KP hierarchy is formulated in terms of the pseudo-differential operator 

L = d + q^ 1 + q 2 d- 2 + q 3 d- 3 + ■ ■ ■ . (5.1) 
The KP flows are defined to be 

d tr L = , L], (5.2) 

with r = 1, 2, 3, ••• (see for instance [35]). These equations admit a Hamiltonian 
description with 

H r = - resU (5.3) 
r 

being the Hamiltonian density of the r-th flow. We have explicitly 

U x = q u (5.4) 

H 2 = q 2 + \ dqi, (5.5) 

n 3 = q 3 + qi 2 + d(q 2 + ^d qi ), (5.6) 

n 4 = q 4 + 3 qi q 2 + - A 8{2q 3 + q x 2 + dq 2 + \ d 2 qi ) (5.7) 
4 o 

and so on. Moreover, when the flows (5.2) are written in terms of the composite fields 

7i, they become 

d tr n 3 = de rd , (5.8) 
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where rj - are composite fields of {q,i\ depending on r, j. It is for this reason that the 
charges Hj = J Hj with j > 1 are all in involution and the KP system is integrable. 

The formulation (5.8) of the KP flows should be compared with the local conservation 
laws of the perturbed parafermion theory. In the parafermion case, h!p — h!p plays the 
role of Tij and d replaces dt r . Therefore, it is natural to identify h!p — hf^ with the 
Hamiltonian densities of the KP hierarchy. Of course, the equivalence between the two 
theories is not exact, because in the KP case the G's can be rewritten as local functionals 
of the H's, while in the parafermion model this is not so. To describe the KP-like 
structure of the parafermion currents, we identify the generators W s as 

W s = , s = 2, 3, 4, • • • (5.9) 

and postulate 

hf - hf = -L_W,, (5.10) 

This correspondence is understood modulo total derivative terms, which are irrelevant 
anyway in both the KP flows and the parafermion conservation laws, since they lead to 
trivial field redefinitions. 

To illustrate the validity of eq.(5.10), we have to use the right normalization for the 
generators W s . Following earlier work on the symmetry [29-31], we scale W s , as 
given by eq.(4.37), by multiplication with 

and choose q = —1/4. Then, in terms of the rescaled variables, eqs.(4.41)-(4.44) become 

hf - hf = W 2 , (5.12) 



* 2 -,o 2 --W 3 , (5.13) 

- # = \ (W 4 + W 2 2 + A d *w 3 ) : (5.14) 

- hf = I (W 5 + 3 W 2 W 3 + - d 2 W 3 ) (5.15) 
8 7 

and so on. Comparison with the expressions (5.4)-(5.7) shows that up to total derivative 
terms (which are gauge dependent in the abelianization method, anyway), eq.(5.10) is 
correct provided that W s = q s -i after the rescaling. We have verified eq.(5.10) for higher 
values of j as well. 

The relation we find between the currents of the perturbed parafermion theory and 
the Hamiltonian densities of the KP hierarchy, generalizes the results we already know 
for the conservation laws of the sine-Gordon model. The sine-Gordon model arises as a 
special case of the perturbed parafermion theory, obtained for u — u. As we have already 
mentioned, this model describes only the neutral sector of our general theory, since both 
components J, J of the U(l) current are identically zero for u — u. All higher spin 
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currents W s with odd values of s also vanish in this sector, as it can be readily seen from 
eq.(4.37). The only non-trivial conservation laws have even spin. Setting u = u = cos8, 
we obtain ip + = ip_ = —06. The conservation law (4.34) for j = 2 reduces to 

d ({d6) 2 ) + g d(cos 2 6) = 0. (5.16) 

For j = 4, eq.(4.46) becomes 

d (5(def - 3{d 2 ef + 2 ded 3 e) + g d (3 cos26(def - sinie d 2 e) = 0. (5.17) 

Similar expressions result for (even) higher values of j, which are all consistent with the 
sine-Gordon equation (1.7). 

It is well known that the currents h^p — of the sine-Gordon model, coincide with 
the Hamiltonian densities of the KdV hierarchy (see for instance [8]). In this case, the 
pseudo-differential operator L in eq.(5.1) satisfies the condition L 2 = d 2 + v and g« are 
all functionally dependent on the KdV field v. We have explicitly q\ = v/2, q 2 = —dv/4, 
q 3 = (d 2 v — v 2 )/8, etc. The relation between v and 9 is given by the Miura map. Flows 
with odd values of j are trivial, since Hj become total <9-derivatives. For the even 
values of j, the KdV reduction of eq.(5.10) describes in fact the relation between the 
conservation laws of the two integrable theories. 

The validity of eq.(5.10), not only in the neutral sector u — u, but in the full per- 
turbed parafermion theory with arbitrary U{\) charge, suggests that there should be a 
2-component reduction of the KP hierarchy, which actually describes the structure of the 
currents fjp — completely. The 2-component reduction in question is the non-linear 
Schrodinger hierarchy, as it has been pointed out recently [36], using a somewhat differ- 
ent method. The Lagrangian description of the parafermion coset perturbed by the first 
thermal operator is the essential ingredient that links our work with theirs. The non- 
linear Schrodinger hierarchy is practically the same as the 2-boson realization of the KP 
hierarchy [37]. The latter was defined following earlier work on the bosonic realization 
of IVoo-type algebras [38] . 

We also point out for completeness that a close relation between the non-linear 
Schrodinger equation and vortex dynamics was discovered several years ago [39], along 
different lines; vortex filaments were intepreted as solitons. This result should not be 
surprizing, given the fact that the perturbed parafermion coset (with g > 0) in 2-dim 
Minkowski space describes the relativistic motion of vortices in constant external field, 
via the Lund-Regge formalism. We think, however, that the role of the non-linear 
Schrodinger equation in vortex dynamics should be investigated further and in connec- 
tion with perturbed conformal field theories. 

We note finally, as a side remark, that it would be interesting to examine the Hamil- 
tonian structure of the local conservation laws of the perturbed coset model. In analogy 
with the KP hierarchy, there might exist a Poisson bracket so that dAj + dBj = could 
be written as 

8A j = {Sl j ,A j }, (5.18) 
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when j > 1 and g ^ 0, for appropriately chosen Hamiltonian Qj. Qj might have a 
natural interpretation in the SU(2)/U(1) coset model, which could clarify the structure 
and physical interpretation of the higher spin conservation laws even further. 

6 Discussion and Conclusions 

In this paper we have studied systematically the Lagrangian description, physical inter- 
pretation and geometry of the parafermion coset model perturbed by its first thermal 
operator. The present work is entirely confined at the classical level and in the large 
iV limit of the theory. Quantum mechanical aspects and 1/N corrections are interesting 
to consider and we hope to return to them in a separate publication. The main tool 
in our study has been provided by the geometric framework of Lund and Regge, upon 
analytic continuation of the coset theory in 2-dim Minkowski (base) space. Reviving 
these geometric ideas in present day conformal field theories seems to be advantageous 
for understanding the Lagrangian description of their integrable perturbations. General- 
izations of these results to other coset models are certainly of great interest. It seems that 
the main problem in the general case is to understand the zero curvature formulation of 
the gauged WZW model (with or without perturbations) and its geometric interpretation 
as Gauss-Codazzi integrability conditions for an embedding problem. 

Many of the issues considered here have been studied before, but with different moti- 
vations in mind. Their relevance in 2-dim field theory and string theory might be quite 
general, not limited to the present models. We expect that the correspondence between 
different string backgrounds and integrable 2-dim field theories could be developed fur- 
ther. Then, a systematic approach to the problem of hidden symmetries in string theory 
might result, as advocated in ref.[15]. We conclude with some applications of the results 
we have obtained so far. 

The infinitely many conservation laws of the perturbed parafermion model have a 
natural geometric interpretation in terms of embeddings. To treat all cases together, 
irrespectively of the sign of the coupling constant g, we summarize the expressions (3.12) 
and (3.27) for the extrinsic curvature as follows 

K + = 2 (K aa + K aT ) = % = - 2kc I u | VI - | u | 2 , (6.1) 

I u I \/l — I u I 



K_ = 2 {K aa - K ar ) = % — _ 2 -2kc\u\ Vl-M 2 - (6-2) 

| u | y 1 — | u | 

The parameter k is a step function: for g < it is k = 0, while for g > we have 
k = 1 and c 2 = g. The metric (3.1) of the 2-dim Euclidean surface S that arises in the 
geometric description of the theory via embeddings, has components 

g aa = | u | 2 , g TT = 1 \ u | 2 . (6.3) 



23 



It is also convenient to introduce the notation g_ = g aa — g TT . 

It is easy to verify that the conservation law of the U(l) current (4.29) reads as 

5 ^ K++Kcg _y g ^ K _ +Kcg _y . (6 . 4) 

The conservation law of the energy-momentum tensor becomes 

d(^K + + 2/tcyW) 2 + l - + 2<? dg- = 0, (6.5) 

while similar (but considerably more complicated) expressions can be obtained for the 
higher spin currents by rewriting the generators W s in terms of the extrinsic curvature and 
the metric. The calculation is straightforward and yields step by step the infinitely many 
integrals of the underlying embedding problem. Their dependence on the original field 
variables is obtained using the defining relations for the metric and the extrinsic curvature 
tensors in terms of r)i(a, r) or Xi(a, r), depending on the sign of the coupling constant g. 
This geometric interpretation is advantageous for understanding the physical meaning of 
the infinitely many charges associated with the classical evolution of the reduced 0(4) 
non-linear <7-model or the string propagation in a 4-dim axionic background of vortex 
type respectively. 

For g = 0, the currents W s become chiral and generate the W^ algebra. This acts 
as a hidden symmetry on the transverse modes of the 4-dim free bosonic string in the 
orthonormal gauge. More details on W^ Backlund transformations in string theory can 
be found in ref. [15]. The lesson for string theory (with or without axionic background) 
is that the classical description of its transverse modes as a perturbed SU{2)/U{\) coset 
model (with g > or g = respectively), implies an infinite number of conservation laws 
at the classical level. It would be interesting to study the quantum mechanical implica- 
tions of this result and try to generalize the correspondence between integrable systems 
living in the string world-sheet and string dynamics, to more arbitrary backgrounds in 
four and higher dimensions. Also, the relevance of this formalism to the quantum theory 
of vortices in (relativistic) superfluids should be investigated further, using the perturbed 
gauged WZW model. Generalization to other integrable perturbations and higher rank 
coset models should be interesting as well. 

Finally, it is straightforward to apply our results to the reduced U{2) Gross-Neveu 
model. As was explained in the introduction, in the special frame (1.8), (1.9), the local 
conservation laws of this model can be obtained from the perturbed parafermion coset 
by the simple substitution (1.10). Generalization to the U(N) case is an interesting 
problem which might be related to integrable perturbations of Grassmannian coset mod- 
els, SU(2)/U(1) being the simplest one. The 1/N expansion of the Gross-Neveu model 
is particularly interesting in this framework, since %p a %l) a develops a non-vanishing vac- 
uum expectation value which breaks chiral invariance. This result might have a useful 
interpretation in quantum coset model calculations. 
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APPENDIX 



In this appendix we summarize the results of the iteration of the recursive relations 
(4.25)-(4.27) for the first few h!p — hf\ in terms of the matrix elements of A . 

- h f = AfA^Af + \ (A] 2 dAf - Af dAl% (A.l) 

fc?> - h ? = \ A^Af (A^Af - A{Aff) - I {A™ d>Af + Af d^) 

- l -Af{A^dAf-AfdA^\ (A.2) 

fl) _ ,(2) _ ,11 ,12.21 (( A\-\\ 2 _ 3 *12 A 2l\ , j* 4 12 ,21/ ,21 o ,12 _ ,12 o.21x 



- fc^ = AqMq ( (Air - - A AqAqJ + ^ 9< - dA£) 



+ \ Af (A? d(AfAf) - Af W)) + 1 d 3 < - Ag 1 d 3 Af) 



+ 1 (2« d 2 < + 3(dAf)d(A^Af) + 3Af(Af d 2 Af + < 9 2 <))(A.3) 
Similarly for tj we have 

4 2 = I AfA? - \ dA™, (A.4) 
tf = ~2 ^^^o 1 ~~ ^ 9Af, (A. 5) 

4 2 = ~^ ((^o 1 ) 2 - ^o 2 ^ 2 ) 1 ) + ^Af + ^Af ~ ^ 2 ^o 2 > (A.6) 
tf = \Af [(Aff - i«) + 1^ fctf + l -Af dAf + ^Af, (A.7) 



tf = \ AfAf ((^) 2 - \a*J%) - \ Af d(AfAf) 

+ l(3Afd 2 Af + 3(dAf)(dAf) + Afd 2 Af), (A.8) 



tf = ~\ AfAf ((Af) 2 - ^4?) - \ Af 9(^M 21 ) 
+ ^ ^ 2 + 4^ 2 &4 21 ) - 1 d 3 A 21 

" \ K ^ + 3«) {dAf) + 3^ 9 2 <) . (A.9) 
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[B^h] + [B ,t 2 ]+t l [t 1 ,B ]-dt 2 + t l dt 1 , (A.10) 

= [B , t 3 ] + [Si, t 2 ] + ti ([t 2 , S ] + [h, S x ]) 

dts + hdh + ^-it^idt. + ^Bo]), (A.ll) 

[S , U] + [Si, i 3 ] + h [t 3 , S ] + *i [t 2 , Si] + ti^ 3 

(* 2 - (^) 2 ) ([t 2 , b ] + Si] + a* 2 ) - du 

(t 3 - ht 2 - t 2 h + (tx) 3 ) ([tx, S ] + dh). (A.12) 
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